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Abstract 

We develop an analytical model of correlated two-electron ionization in strong infrared laser 
fields. The model includes all relevant interactions between the electrons, the laser field, and the 
ionic core nonperturbatively. We focus on the deeply quantum regime, where the energy of the 
active electron driven by the laser field is insufficient to collisionally ionize the parent ion, and the 
assistance of the laser field is required to create a doubly charged ion. In this regime, the electron- 
electron and the electron-ion interactions leave distinct footprints in the correlated two-electron 
spectra, recording the mutual dynamics of the escaping electrons. 
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I. INTRODUCTION 



In strong infrared laser fields, following one-electron ionization of an atom or a molecule, 
the liberated electron can recoUide with the parent ion [1, 2]. The electron acts as an "atomic 
antenna" [2], absorbing the energy from the laser field between ionization and recoUision 
and depositing it into the parent ion. Inelastic scattering on the parent ion results in 
further collisional excitation and/or ionization. Liberation of the second electron during the 
recoUision - the laser-induced e-2e process - is known as correlated, or nonsequential, double 
ionization (NSDI). 

The phenomenon of NSDI was experimentally discovered by Suran and Zapesochny [3] 
for alkaline-earth atoms (for further experimental investigations of NSDI for alkaline-earth 
atoms, see, e.g., Refs. [4-6]). In this case, autoionizing double excitations below the second 
ionization were shown to be extremely important. For a theoretical study of these effects, 
see, e.g., Ref. [7]. For noble gas atoms, nonsequential double ionization was first observed 
by L'HuiUier et al. (see, e.g., Refs. [8, 9]). The interest to the phenomenon of NSDI 
grew rapidly after NSDI was rediscovered in 1993-1994 [10, 11], and now for IR fields and 
higher intensities. Recently, correlated multiple ionization has also been observed [12, 13]. 
The renewed interest in NSDI has been enhanced by the availability of new experimental 
techniques that allow one to perform accurate measurements of the angle- and energy- 
resolved spectra of the photoelectrons, in coincidence. Such measurements play a crucial 
role in elucidating the physical mechanisms behind the NSDI. 

Prom the theoretical perspective, direct ab initio simulations of the photoelectron spectra 
corresponding to correlated (or nonsequential) double ionization in intense low-frequency 
laser fields represent a major challenge. Only now such benchmark simulations have become 
possible [14] for the typical experimental conditions (the helium atom, laser intensity / ~ 
10^^ W/cm^, laser wavelength A = 800 nm). 

What are the physical reasons behind these numerical challenges, which push the modern 
computational resources to their limit, occupying thousands of processors for weeks at a 
time? First, they lie in the need to deal with highly nonstationary two-electron dynamics, 
with electron energies changing by hundreds of eV on a subfemtosecond time scale, the 
characteristic amplitudes of electron oscillations reaching several tens of angstroms, and 
final electron energies ranging from zero to 10^ eV. Accurate description of such dynamics 



3 



requires attosecond-scale time steps, very large grids, and small grid steps ~ 0.1 A. 

Second, one needs to analyze the results of such massive calculations, which output a five- 
dimensional, time-dependent wave function (one spatial dimension is saved by the cylindrical 
symmetry of the problem in a linearly polarized laser field.) Extracting essential physical 
processes and mechanisms responsible for correlated double ionization from such massive 
data arrays is a separate and equally formidable challenge. 

In addition to the experimental measurements, especially of the correlated electron spec- 
tra [12, 13, 15-19], tremendous insight into the physics of the problem has been obtained 
from classical simulations performed in Refs. [20-25]. These papers have demonstrated 
a variety of the regimes of nonsequential double and triple ionization. Not only do these 
simulations reproduce key features observed in the experiment, they also give a clear view 
of the (classical) interplay between the two electrons, the potentials of the laser field, and 
of the ionic core. They also show how different types of the correlated motion of the two 
electrons contribute to different parts of the correlated two-electron spectra. 

Our goal is to develop a fully quantum, analjd;ical treatment of this problem. It is 
important that the approach takes into account all relevant interactions - those with the 
laser field, the ion, and between the electrons - nonperturbatively. Given the complexity 
of the problem, it is clear that the analytical description will have to incorporate physical 
understanding of the dynamics gained from the previous experimental and theoretical work. 

In particular, the physics of double ionization is different for different intensity regimes, 
separated by the ratio of the energy of the recoUiding electron to the binding (or excitation) 
energy of the second electron, bound in the ion. 

The maximum energy, which the recoUiding electron can acquire from the laser field, is 
~ 3.2Up [1], where Up — F'^/Acu'^, F is the laser field strength, and uj is the laser frequency 
(atomic units are used throughout the paper). Even when 3.2C/p is far from sufficient to 
hberate other electrons, experiments have observed correlated ionization [12, 13, 17-19]. 
As opposed to the more conventional high-C/p regime (see, e.g., Refs. [15, 16, 26-35] and 
references therein), in the low-Up regime the assistance of the laser field during the recoUision 
is crucial. 

Here, we focus on this most challenging \ow-Up regime, where the nonperturbative in- 
terplay of all three interactions is crucial. In this regime, existing classical and quantum 
analysis (see, e.g., Refs. [23, 25, 36]) demonstrates two possibilities of electron ejection after 



4 



the recoUision. First, the two electrons can be ejected with httle time delay compared to the 
quarter-cycle of the driving field. Second, the time delay between the ejection of the first 
and the second electron can approach or exceed the quarter-cycle of the driving field. In 
these two cases, the electrons appear in different quadrants of the correlated spectrum. If, 
following the recoUision, the electrons are ejected nearly simultaneously, their parallel mo- 
menta have equal signs, and both electrons are driven by the laser field in the same direction 
toward the detector. If, following the recoUision, the electrons are ejected with a substantial 
delay (quarter-cycle or more), they end up going in the opposite directions. Thus, these two 
types of dynamics leave distinctly different traces in the correlated spectra. 

Here, we consider the case in which the two electrons are ejected nearly simultaneously. 
We show that in this case the correlated spectra should bear clear signatures of the electron- 
electron and electron-ion interactions after ionization, including the interplay of these inter- 
actions. We identify these signatures. 

To study NSDI analytically, we use the strong-field eikonal-Volkov approach (SF-EVA) 
and follow the recipe described in Ref. [37]. Our model complements earher theoretical 
work on calculating correlated two-electron spectra using the strong-field S-matrix approach 
[26] . The key theoretical advance of this work is the ability to include nonperturbatively all 
relevant interactions for both active electrons: with each other, with the ion, and with the 
laser field. Electron-electron and electron-ion interactions are included on an equal footing. 
Our model ignores multiple recoUisions and multiple excitations developing over several laser 
cycles, such as those seen in the classical simulations [22]. This simplification is particularly 
adequate for the few-cycle laser pulses, as demonstrated in the experiment [38] . According 
to this experiment, multiple recollisions are noticeably suppressed already for 12 fs pulses 
at A = 800 nm. For 6-7 fsec pulses at A = 800 nm, such simplification is justified. 

The rest of this paper is organized as follows: In Sec. II, preliminarily comments from the 
point of view of the Landau-Dykhne adiabatic theory are made. Basing on the SF-EVA, we 
develop the analytical model of NSDI in Sec. III. The correlated spectra are first calculated 
within the strong field approximation (SFA) in Sec. IV. In Sees. V and VI, the roles of 
electron-electron and electron-ion interactions arc analyzed. We show how these interactions 
lead to profound qualitative differences from the SFA-based models in positions and shapes 
of maxima in the correlated spectra. Finally, the conclusions are made in the final section. 
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FIG. 1: The diagram of NSDI within the considered regime. 



II. FUNDAMENTAL PRINCIPLES AND BACKGROUND 

The process is shown by the Feynman diagram in Fig. 1. The system begins in the 
ground state \gg) at time tj. At an instant t;,, intense laser field promotes the first electron 
to the continuum state |k); the second electron remains in the ground state of the ion l^'"''). 
RecoUison at frees both electrons. The symmetric diagram where electrons 1 and 2 change 
their roles is not shown, but is included in the calculated spectrum. 

Despite the fact that the main approach used in this paper is the strong-field eikonal 
Volkov approximation (SF-EVA), it is methodologically useful to look first at the problem 
at hand from the point of view of the Landau-Dykhne (LD) adiabatic approximation. Indeed, 
since the frequency of the laser field is low compared to other relevant energy scales such 
as the binding energy, the LD approximation is a natural way of tackling the problem. 
Moreover, the LD approach is an especially handy tool for obtaining results within the 
strong field approximation, which is used as the zero-order approximation in the SF-EVA 



According to the LD method [39, 40] (see also Ref. [41]), if the Hamiltonian of a system 



probability F of the transition ipi ^ ipf is given by (the atomic units are used throughout) 



[37]. 



H(t) is a slowly varying function of time t, and H{t)ipn{t) = En(t)ipn(t) {n = i, /), then the 




(1) 



where ti is any point on the real axis of t, and to is the complex root of the equation 



(2) 



which lies in the upper half-plane. If there are several roots, we must choose the one that is 
the closest to the real axis of It must be stressed that there are no assumptions regarding 
the form of the Hamiltonian. Further discussions and generalizations of the LD method can 
be found in Ref. [42-46]. 

The LD approach has many applications in different areas of physics. In particular, it 
has been extensively used in strong field physics [41]. For example, the problem of single- 
electron ionization can be analyzed within the LD approximation (see Ref. [41, 47-51] 
and references therein) by setting Ei — —Ip to be the energy of the ground state {Ip is 
the ionization potential) and Ef{t) — [k -|- A(t)]^ /2 to be the energy of the free electron, 
oscillating in the laser field. Here k is the final momentum of the electron at the detector, 
and 



is the vector potential. The results of the LD approach in this form are fully consistent with 
the usual SFA. Improving the SFA result by incorporating the Coulomb potential is also 
possible, as we discuss below for the case of two liberated electrons. 

Let us apply the general approach Eq. (1) to the two-electron process under consideration. 
The NSDl has two stages, namely ionization of the first electron and the rccoUision. Hence, 
strictly speaking, the LD adiabatic approximation has to be applied to each of the two stages, 
since the total amplitude of the process is the product of the ionization amplitude and the 
recoUision amplitude. However, it is the second (recoUisioin) amplitude that is responsible 
for the shape of the correlated spectra. The first amplitude only gives the overall height of 
the spectra, as it determines the overall probability of the recoUision. Since at this stage we 
are only interested in the shape of the correlated spectra, we omit the ionization amplitude 
from this discussion (it is included later in the full treatment). 

As a zero approximation, we define Ei{t) and Ef{t) for the second part of NSDI without 
the Coulomb interaction. Before the recoUision at the moment tr (Fig. 1), one electron is 
bound and another is free. The classical energy of the system before tr is 



A{t) 



(F/lu) sm{ujt) 



E,{t)^^[-A{t,{t)) + A{t)f + E,+ , 



(3) 
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where Eg+ denotes the energy level of the second (bound) electron. The time of "birth" 
(ionization) for the first electron tb{t) is the standard function of the instant of recoUision 
tr, which is obtained from the saddle-point S-matrix calculations in the Appendix. In Eq. 

(3), we have assumed that the recolliding electron has been born at tiy{tr) with zero velocity. 
After the recoUision, both electrons are free and the energy of the system is 

Ef{t) = \[k,+A{t)r+l[k,+A{t)r, (4) 

where ki 2 are the asymptotic kinetic momenta at i ^ 00 of the first and second electrons, 
corr esp ondingly. 

Now, substituting Eqs. (3) and (4) into Eq. (1), we obtain the correlated spectrum 
standard for the strong field approximation (SFA), 

r5F^(ki,k2) pa exp ( --Im5'5F^(ki,k2) ] , (5) 



UJ 



SsFA{ki,k2) 



,- 

\ (ki + A{^)f + ^ (k2 + A{^)f - \ [A(^) - A($(7; ^))f + /f 



where the phase of "birth" (ionization) $(7; v?) corresponding to the rccollison phase if and 
the transition point are defined by Eqs. (All) and (A8) in the Appendix. We will come 
back to more rigorous analysis of the same spectra in the next section. 

The major stumbhng block is to account for the electron-electron and the electron-ion 
interactions on the same footing, nonperturbatively. To include these crucial corrections, 
we have to include the corresponding Coulomb interactions into Eij(t). With the nucleus 
located at the origin, the electron-electron and the electron-core interaction energies are 

Ke = l/|ri2(t)|, Vj^''^ = -2/\r,,,{t)\, (6) 

correspondingly. Here ri2(t) = ri{t) — r2{t) and ri_2(i) are the trajectories of the two 
electrons. 

However, we immediately see problems. The corrections depend on the specific trajectory, 
and one needs to somehow decide what this trajectory should be. Note that the classical 
trajectories ri 2(t) in the presence of the laser field and the Coulomb field of the nucleus may 
even be chaotic. The solution to this problem has already been discussed in the original 
papers by Popov and co-workers [52-55] for single-electron ionization. In the spirit of the 
eikonal approximation, these trajectories can be taken in the laser field only [37, 54, 55], so 



that they correspond to the saddle points of the standard SFA analysis. Not surprisingly, 
in the SFA these trajectories start at the origin, 

J to 

However, here we run into the second problem: the potentials Vge and Ven' are singular. 
Consequently, the integral in Eq. (1) is divergent and the result is unphysical. Therefore, 
such implementation of the Coulomb corrections requires additional care. 

The next sections describe a rigorous approach that deals with these two problems, both 
defining the relevant trajectories and removing the divergences of the integrals. 



III. BASIC FORMALISM 



The key step in dealing with the singularities of the Coulomb potentials during the rec- 
oUision is to partition the electron-electron and electron-ion interactions in the two-electron 
Hamiltonian as follows: 

^ee — ^ee ^ee,lng ~l~ ^ee,lng ^ee,lng ~l~ ^^e,s?jr) 

^en = ^en ^en,lng ~l~ ^en,lng ^en,lng ~l~ ^^en,shr- (7) 

The potential Vexing has a long-range behavior identical to T4e, but no singularity at the 
origin, and l^Vee,shr is singular but short-range potential. The same applies to Ven,ing and 

^Ven,shr- We choOSC 

^Ven,shrir) = Ven{r)exp{-r/ro), AVee,shriri2) = Vee{ri2) -ru/r^S , (8) 

where tq and rf^2 ^i^l be defined later. Note that the partitioning (7) and (8) has been 
employed originally in the Perelomov- Popov- Terent'ev approach [52-55] for the problem of 
single-electron ionization. 

Now, we can write the Hamiltonian as 

H{t) = Hs{t) + AVshr, 

where AVshr = AVee,shr + AVen,shi- and Hs is the rest, which includes smoothed Coulomb 
potentials for electron-electron and electron-nuclear interactions. Vexing and Ven,ing- 
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To first order in AVghr, the amplitude to find two electrons with momenta ki,k2 at the 
detector at the time t is 

a(ki, ka) = -i dtr j (fk (kiksl Usit, tr)AVshr |^+k> (k^+| C/,(i„ U) \gg) . (9) 

Approximations in Eq. (9) are first order in AVshr and the assumption that at the moment 
of recoUision the ion is in its ground state. Both are well justified. 

The next step is to approximate the two parts of the evolution: before tr and after tr- 
The key component for correlated spectra is the second part - after tr- The main aspect of 
the first part of the evolution - prior to U - is to supply an active electron with the required 
energy. 

To simplify the amplitude {]s.i]s.2\ Us{t,tr) AVghr \g~^k) , we insert the decomposition of 
unity, 

b{ki,k2.k,tr) = {kik2\Us{t,tr) AVshr \g^k) = 

^ J j dhidh2 (kikal Us{t, tr) |rir2) (riral A^^^ |^+k> ^ 

^ J J dhidh2 (ki + A{tr), k2 + A{tr)\ Tirs) (riral AV^hr |^+k> x 

exp [-^ {i[ki + A(r)]2 + l[k, + A(r)]2+ 

+Vee,ln9i^l2{T)) + Ven,lngi^li^)) + VenMgi^^))] dr ■ (10) 

Here we have applied the SF-EVA method [37]. The integral from the nonsingular parts 
of the electron-electron and electron-ion interactions are calculated along the trajectories in 
the laser field only. The trajectories 

riAt) = ri,2 + / [ki,2 + A{T)]dr (11) 

and ri2{t) — ri{t) — r2{t) begin at the positions ri, r2 at instant tr- The bra- vectors 
(ki^2 + A(ir) I are plane waves. Their distortion by the electron-electron and electron-core 
interactions appears in the (ri, r2) -dependent exponential phase factors in Eq. (10). 

Since AVghr is a short-range potential and l^'"'") is limited within a characteristic ionic 
radius, the term (rir2| AVghr \9~^k) allows us to fix the initial values of ri and r2. The 
characteristic radius for the partitioning of the Coulomb potential into the short-range and 
long-range parts is set as ro = r^2 —'^1 1 -^5+ 1- Therefore, we pull the exponential factor out 
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of the integral in Eq. (10) with ro = r^"^ — 1/ \Eg+\ and r\ — r2 — 0, 

-^y^*|i[ki + A(r)f+ 

(12) 



6(ki, ka, k, U) ^ (ki + A{tr) ka + A{tr)\ AVshr \g'^k.) exp 

+ ^[k2 + A(t)]2 + VeeMgiT^uir)) + VenMgiM^)) + Kn.Zng (r2 (t) ) | rfr 

Effects of the long-range tails of Vee and Ven appear in the exponent while the coUisional 
transition is govered by the short-range interaction [37]. The states |k), so far, represent 
any convenient basis set of continuum sates in the laser field. 
To simplify the amplitude 

c{'k,tr) = {kg+\ Us{tr,ti) \gg) , 

we note that the second electron is bound during the whole evolution, and hence we can 
simplify c(k, ir) using single active electron approximation. In this approximation, Us{tr,ti) 
describes one-electron dynamics in the self-consistent potential of the ionic core. 

The effective Hamiltonian for evolution between ti and tr is 

^,e(ri, t)=ki + Kc(ri) + VL{ri,t), 

where Ki is the kinetic energy operator and VL{ri,t) is the interaction with the laser field. 
Now the amplitude c(k, tr) becomes 

dtb (k| Usc{tr: 4) VL(ri, tb) Igo) exp [i \Egg\ {t^ - t^)] , (13) 

where Iqd) = {gt \gg) is proportional to the Dyson orbital between the ground states of the 
neutral and ion. 

The ionic potential contributes to the propagator in Eq. (13) twice: when the electron 
leaves the atom near tb and when it returns to the ionic core near tr- The contribution "on 
the way out" introduces standard Coulomb correction [52-57] to the ionization amplitude 
and hence affects the overall height of the final two-electron distribution. The contribution 
of Vsc "on the way in" affects the spatial structure of the recoUiding wave packet. As shown 
in Refs. [37, 58], for short collision times the Coulomb-laser coupling is small and Vgc "on 
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the way in" can be included in the adiabatic approximation, 



C{k,tr) \Kvg^) ~ -iRc / dtb \Kvg^) {k+A{tb) - A{tr)\ VL{tb) \9d) X 



exp 



[k+ A(r) - A{tr)fdT + i \Eg+\ {tr -tb)+i \Egg\ {tb - ti) 



(14) 



Here |ke^,) is the field-free continuum wave function in the eikonal approximation, 
which includes distortions of the incoming plane wave with asymptotic momentum k, 
(k-|- A(i;,) — A{tr)\ is a plane wave, and Rc is the Coulomb correction to the ionization 
amplitude which compensates for approximating (k-|- A{tb) — A{tr)\ with a plane wave in 
the matrix element (k + A{tb) — A{tr) \ 14 (t^) \gD)- 

Now, putting together Eq. (14) and Eq. (12) and changing the integration variable 
k — >• k + A(tr), we arrive at 

a(ki,k2) ^ ~ j dtr J dtb J d^k J dhid^r2 

12, 



X exp 



[k + A{T)YdT + i \Eg+ 1 {tr -tb) + i \Egg\ {tb - ti)- 



tb 



-t / I ^[ki + A(r)]2 + ^[k2 + A(r)]2 + Ke ,ing(ri2 ("?■)) + K:n,ing(ri(T)) + Ven,lng{t^2 (r))|cir 
X (ki + A{tr), k2 + A{tr) I rira) (rir2| AVshr |^+, ke^ + A(t,)> Rc (k + A(4)| ^^(4) Igo) • (15) 

Note that if the Coulomb corrections Vee,ing and Ven,ing ai'e ignored in the exponent of Eq. 
(15), then Eq. (15) coincides with Eq. (Al) within exponential accuracy. 



IV. THE CORRELATED TWO-ELECTRON IONIZATION WITHIN THE SFA 

In this section, we find the correlated spectrum of the NSDI by using the strong field 
approximation (SFA). In the next sections, we will improve the SFA result by employing the 
SF-EVA [37], i.e., the perturbation theory in action with the SFA result as the zero-order 
approximation. 

Ignoring the Coulomb corrections in Eq. (15) and performing the saddle-point calcula- 
tions described in the Appendix, we reach the usual SFA expression for the correlated NSDI 
spectra - Eq. (5). 

To illustrate the SFA results, we plot the two-electron correlated spectrum for a system 
with Ip of Ar in Fig. 2. In Fig. 2.A we set 7 = 0. Such an SFA spectrum has a peak at 
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FIG. 2: Correlated two-electron spectra (5) of Ar (linear scale) within the SFA at 7 x 10^^ W/cm^, 
800 nm = k^2 — 0) {^) 7 = 0; (b) 7 = 1.373. Maxima of figures correspond to probability 

densities: (a) 1.7 x 10"^ (b) 2.9 x lO'^^. 



^Ili ~ ^Ib ~ "~0.78 a.u., which is the maximum of the vector potential —F/uj ~ —0.78 a.u. 
The last fact has the following interpretation: NSDI is most efficient when the velocity of the 
incident electron is maximal. This is achieved near the zero of the laser field, = F cos if, 
and the maximum of A((y9). An electron liberated at this time could acquire the final drift 
velocity ~ —F/uo. However, including the correct value of the Keldysh parameter 7 not 
only substantially shifts the peak position (Fig. 2.B), but also lowers the maximum by nine 
orders of magnitude. 



V. ELECTRON-ELECTRON INTERACTION 



In this section, we demonstrate the changes in the correlated spectrum due to the electron- 
electron repulsion. 

Coulomb corrections to the single-electron SFA theory were first introduced by Perelomov, 
Popov, and Terent'ev [54, 55] using the quasiclassical (imaginary time) method (for reviews, 
see Refs. [59, 60]). More recently, further improvements to this method have been considered 
in Refs. [56, 57]. These improvements considered not only subbarrier motion in imaginary 
time, but also the effects of the Coulomb potential on the phase of the outgoing wave packet 
in the classically allowed region. These improvements allowed the authors of Refs. [56, 57] 
to obtain quantitatively accurate results not only for ionization yields, but also for the 
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above threshold ionization spectra of direct electrons (i.e., not including recoUision) . An 
alternative, but conceptually similar, approach is the SF-EVA [37] . Unhke the two previous 
methods, the SF-EVA allows a simple treatment of the electron-electron and electron-ion 

interaction in the two-electron continuum states. 

According to the SF-EVA, the contribution of the interaction potentials is calculated 
along the SFA trajectories. 



ri,2(<^) 



1 



[ki,2 + A(0)]d0. 



Note that at the moment of recollision ip^, the electrons are assumed to be at the origin, 
ri,2(<^r) = 0. However, this does not cause any divergence since according to Eq. (15) we 
have to use the regularized potential Vee,ing- 

From Eqs. (7) and (8), the potential energy of electron-electron repulsion along these 
trajectories is given by 



1 — exp 



ri2{^) = \ri{ip) - r2{(p)\ 



'12 



H 2 



+ 



{k±i - k±2) 



UJ 



(16) 



As discused in Sec. Ill, the parameter is set to the ionic radius, — 



The correlated spectrum, which accounts for the electro-electron interaction, has the form 



ree(ki,k2) pa exp 



5'ee(ki,k2) 



Re<^o 



--Im [SsfaO^I: ka) + 5'ee(ki, ks) 

CO 



Vee,lng{^)d(p. 



(17) 



Figure 3 shows the contribution of electron-electron repulsion to the spectra of NSDI for 
an atom with Ip of Ar (for experimental data see Ref. [17]). 

Comparing Figs. 2 and 3, we readily notice a dramatic influence of electron-electron 
interaction on the correlated spectra. Electron-electron repulsion splits each SFA peak into 
two peaks because, due to the Coulomb interaction, two electrons cannot occupy the same 
volume. Note that the larger the difference between the perpendicular momenta of both of 
the electrons, the closer is the location of the peaks. 
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VI. ELECTRON-ION INTERACTION 



Now we include the electron-ion attraction. The potential energy of electron-ion interac- 
tion for the case of two electrons and a single core, after partitioning (7) and (8), is 

V^r.,lnM = -2 J] 4- f 1 - exp [-^1 y 




riM = \l{h..^^ + ^2(^os^-cos^^,)y +(^k^,,,^^ (18) 

As far as the parameter Tq is concerned, we set it equal to = V-^p^^- 

Finally, the correlated spectrum of NSDI, which takes into account both electron-electron 
and electron-ion interactions, reads 

ree+en(ki, ka) ^ cxp Im [SsfaO^i, ^2) + 'S'eelki, ka) + k2)]) , (19) 

,0 



5'e„(ki,k2) = / Ven,ln9{'P)d'P 



To illustrate the influence of the electron-ion attraction, wc have plotted the correlated 
spectra of Ar in Fig. 4 for different perpendicular momenta. Comparing Figs. 3 and 
4, we conclude that the larger the difference between the perpendicular momenta of the 
two electrons, the larger is the contribution of the electron-ion interaction. Furthermore, 
accounting for electron-ion attraction increases the probabihty of NSDI by 15 orders of 
magnitude. This occurs because, as in the case of single-electron ionization, electron-core 
interaction signiflcantly lowers an effective potential barrier. We can also conclude that 
correlated spectra pictured in Figs. 4.C, 4.D, and 4.E have the biggest contribution to the 
total probabihty of NSDI, which is an integral of the probabihty density over momenta of 
both of the electrons. Note that, on the one hand, the maximum of the probability density 
shown in Fig. 4.F is the largest among those presented in Fig. 4, and on the other hand, 
this maximum is localized in a few pixels; therefore, the integral contribution of Fig. 4.F 
to the total probability is smaller than Fig. 4.E. Additionally, as one would expect, further 
increasing of k± leads to a decrease of probability density. The correlated spectra in Fig. 4 
agree with the experimental data [17] in quadrants one and three. The considered diagram 
(Fig. 1) does not contribute to signals in quadrants two and four. Note that taking into 
account a nonzero value of 7 is vital to achieve agreement with the experimental data. 
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Prom Eqs. (16) and (18), we can notice that if oo and ro — > oo, the Coulomb 

corrections Veering and Vending vanish, and the SFA result is recovered. Therefore, we conclude 
that the radii r^^' and Tq contain the information about the initial position of electrons after 
they emerged in the continuum. Obviously, the intra-electron distance should be on the 
order of an ion radius. 

VII. CONCLUSIONS 

The analytical quantum-mechanical theory of NSDI within the deeply quantum regime, 
when the energy of the active electron driven by the laser field is insufficient to coUisionally 
ionize the parent ion, has been formulated based on the SF-EVA approach. On the whole, 
the presented model agrees with available experimental data [17]. We have defined the 
quantum-mechanical phase of birth of the active electron (All), which accurately accounts 
for tunneling of the recoUiding electron in the regime where both the phases and are 
complex. Moreover, it has been demonstrated that ignoring such a contribution of tunneling 
of the active electron fails to agree with the experimental data. 

Furthermore, our results show that any attempt to interpret NSDI spectra in this regime 
in terms of a simple SFA-based streaking model would lead to wrong conclusions on the 
relative dynamics of the two electrons. 

The contributions of the electron-electron and electron-ion interactions have been an- 
alyzed. Both play an important and distinct role in forming the shape of the correlated 
spectra. 

The presented model is not able to give the correlated spectra in quadrants two and 
four. It is because the considered process, when two electrons detach simultaneously from 
the atom, does not contribute to that area. However, we incline to believe that those 
parts of the correlated spectra are formed due to recoUision and excitation of the ion plus 
subsequent field ionization [61, 62], and it should be noted that this mechanism has been 
also observed in classical simulations [23]. We are planing to develop an analytical model of 
such a mechanism in future papers. 
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APPENDIX A: THE PHASE OF IONIZATION OF THE RECOLLIDING ELEC- 
TRON AS A FUNCTION OF THE PHASE OF RECOLLISION 

Employing the SEA, we write the formula corresponding to the digram of NSDI (Eig. 1) 

/ dh [ dtr I d''\^U{t,tr)—\\^g+){g+\^\VL{h)\gg)x 

expj-^^ [k + A{T)]^dT + z\Eg+\{tr-tb)+z\Egg\{U-ti)Y (Al) 

where U{t,tr) is the evolution operator of the studied system, ri2 is the distance between 
the electrons, VL(tb) is the interaction between the ionized electron and the laser field, and 
Eg+ and Egg are energies of the states \g~^) and \gg), respectively. 

We use the saddle point approximation (SPA) to calculate the integrals over k and tb in 
Eq. (Al). The phase of the integral over k has the following form: 

The saddle point of this integral is given by 

ft 



k* 



r A{r)dr, (A2) 
ii- H Jtb 



with the restriction tj. ^ tb- Note that generally speaking, k* can be complex since tb, as 
will be clarified below, is complex for 7 7^ 0. The phase of the integral over tb in Eq. (Al) 
reads 

-^2(^6) = / ' + A{T)fdT + \Eg+\{tr - tb) + \Egg\{tb " t i) . 

Hence, the saddle point tbit^) is a function of t^ and given as a solution of the following 
equation 

cos(a;t^) — <zo^\ijjtb{trS\ . , , . , 
I + ^v^{utb{tr)\ = ±Z7, (A3) 

where 7 is the Keldysh parameter for the first electron. 



^^^VSS+HS. (A4) 
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Methods of computing the saddle points [Eqs. (A3) and (A7)] have been widely discussed in 
the literature (see, for example, Refs. [31-33] and references therein). In the current paper, 
we use a general and simple approach for identifying correct saddle points between different 
solutions of the saddle-point equations in the complex plane (see Sec. IV). 

It is convenient to introduce the following phases: 0b = utf, and (pr = ^ir- According to 
Eq. (A3) , the saddle point 06 is a complex double- valued function of which can be given by 
06(+7; Vr) and 06(— 7; </7r-), whcrc 7 is the Keldysh parameter (A4). Here, the complex single- 
valued function 0^(7; </?r) is defined as a solution of the following transcendental equation: 

cos^^ -cos 06(7;^^) . , / N ■ /ACTN 
— hsm06(7;(^r) = Z7. (A5) 

No analytical solution of such an equation is available. The special case of the function 

Mr, Vr) for 7 = 0, 

— ^ hsm(/7b((^^) = 0, (A6) 

is very important because of the following two reasons. 

First, the function ipi)[ipr) is a real valued function for real Lpr (and single-valued for 
any complex (pj.)', this allows one to interpret the motion of the first electron in terms of 
classical trajectories. The function (fbifr) is defined on the interval (7r/2,27r] because only 
during that interval can the free electron recoUide with its parent ion. The function (fib{(fir) 
is bounded in the interval ^ fbifr) < 7''/2. Second, the function (pbifr) can be physically 
understood as a tunneling limit (7 <^ 1) of 0^(7; fr) (i-c, the low-frequency limit). 

In terms of the laser phase, the vector potential A{ip) is 

A{ip) = — (F/o;) simp. 

The difference between Eq. (A5) and Eq. (A6), which connect the phase of recoUision (pr 
with the phase of ionization (pbifr) or 0^(7; (pr), is also important for the last step of NSDI 
- the release of the two electrons following the recoUision at (pr- 

Calculating the integral over tr in Eq. (Al) by the SPA, we need to obtain the transition 
point ip^ for negligible 7. It is the solution of the saddle-point equation 

AE{f'^) ^l[k, + A(^°)]% i [k, + A{p'^)Y - \ [A{fl) - Kipbi^m" + = 0' (A7) 
such that 

7r/2 < Re(/?° ^ 27r, 
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where Ip ' — \Eg+ \ is the ionization potential of the second electron. For 7 7^ 0, the equation 
is 



AE{r, ^l) + H^'r)f + ^ [k, + A{^l)f - \ [A((^°) - A(0,(7; (^^))]' + /f = 0, 



where ^6(7; V'^) now depends on 7. Note that Eqs. (A7) and (A8) are basically Eq. (2) 
written in slightly different notations. 

If the solution of Eq. (A7) on the interval (vr/2, 27r] is real, then direct coUisional ionization 
is possible. However, we are interested in the deep quantum regime when the following 
inequality is valid for the second electron: 



By introducing the Keldysh parameter for the second electron 72 = W we can write the 
last inequality as 



Equation (A9) physically means that the returning electron does not have enough energy 
to free the second electron. 

When recoUision energy is not sufficient for coUisional ionization, transition requires help 
from the laser field. Mathematically, the arising integral is similar to those in the Landau- 
Dykhne (LD) theory (see Sec. II). The energy gap AE{(pr) in Eqs. (A7) and (A8) plays 
the role of the transition energy for the LD transition [the term Ei{t) — Ef{t) in Eq. (1)]. 
The peculiarity of AE{'~f; ipr) given by Eq. (A8) is that it need not be real even for real 
(fr, since in the term [A{(pr) — ^{4>b{l', y^r))]^ the phase 0f, is complex. This subtle aspect 
underscores the important difference between using the solutions of Eq. (A5) or Eq. (A6) 
for the phase of ionization 0^. For classical trajectories, where (pbi'fr) [Eq. (A6)] is real 
for real cpr, AE{(pr) [Eq. (A7)] is also real for real cpr- This is the standard assumption 
in the LD theory. When the complex phase of ionization (pbil] <fr) [Eq. (A5)] is used, i.e., 
when "quantum" trajectories for rccollision arc used, AE{'~f; ip^) [Eq. (AS)] need not be real 
for real (pr- In [45], the LD method has been generalized for this case, provided that the 
complex function A£J(7; ip) satisfies the Schwarz refiection principle (recently, this result 
was confirmed and further generalized in Ref. [46]). 

The function AE{(p) (A7) obeys the Schwarz refiection principle, i.e., AE*{(p*) — AE{(p). 
Hence, we can conclude that if </7° is the solution that hes in the lower half-plane, then (y^^)* 



(AS) 



7^') > 3.17Up. 




72 > 1.26. 



(A9) 
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is the solution that hes in the upper half-plane. Furthermore, it can be easily proven that the 
following equation takes place for any function /S.E{ip) which satisfies the Schwarz reflection 
principle and any complex number ip^: 

Im / AE{(p)d(p = -Im / AE{(p)d(p. 

From the previous equation, we can see that the transition points that lie in the lower 
half-plane lead to exponentially large probabilities, which are unphysical. Hereafter, let (p^ 
denote the solution of the equation AE((p^) — 0, which is the closest to the real axis and 
lies in the upper-half plane. 

Before continuing our discussion, let us point out the following simple equalities, 
which follow from Eq. (A5): Re [05(+7; = Re 7; V^r)] and Im [0b(+7; 99^)] = 
— Im 7; (fir)] for real ip^- Furthermore, we obtain 

rb{-r,^*r)=Mr,^r), Im / ' AE{r,^)d^ = -l^ AE{-r,ip)dip, (AlO) 

where £'(7; </?) is defined in Eq. (A8). 

Bearing in mind that formula (Al) must give an exponentially small result (which implies 
that the transition point must be located in the upper-half plane) and taking into account 
Eq. (AlO), we define the phase of ionization in the case of 7 7^ as 

(f>b{-T, Vr) if Im (</7r) > 0, 

Re [06(7; (/p,)] if Im ((^,) = 0, (All) 

0b(+7; V^r) if Im {ifr) < 0. 

Equation (All) is the most consistent definition of the quantum-mechanical phase of 
ionization of the first electron as a function of the phase of return. Generally speaking, 
there was an ambiguity in selecting the value of $(7; (pr) for real (pr- However, we have 
chosen it in such a way due to the following reason. The function Im [$(7; (/j^)] has a jump 
discontinuity on the real axis, but the function Re [$(7; (fir)] has a removable discontinuity 
that can be removed by employing the equality 

Re [0^(7; (fir)] = ^ [(pbir, 'Pr + iO) + M-T, 'Pr - iO)] (for real (fir). 

Furthermore, the function $(7; (^r) obeys the Schwarz refiection principle [$*(7;<^*) — 
$(7; (fir)], and the following equality takes place: 

^{0;(Pr) = fbifr)- 
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$(7; (fir) 



It is essential that according to definition (All), the function $(7; ipr) is real- valued on the 
real axis and thus allows the identical interpretation in terms of the classical trajectories as 
for ipi,{ipr)- Therefore, the definition of </?^ and inequality (A9) are unchanged in the case of 
7 7^0. 
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FIG. 3: Role of electron-electron interaction. Correlated spectra of Ar (linear scale) at 7 x 10^^ 
W/cm2, 800 nm are calculated using Eq. (17) with ryl = 0.985 a.u. (7 = 1.373). Electron-core 
interaction is not included. Spectra are shown for different values of (in a.u.) for both electrons: 
(a) k^^ = k^2 = 0; (b) fcxi = 0,k^2 = 0.2; (c) k±^ = -0.l,k±2 = 0.2; (d) k^^ = -0.2, k^2 = 0.3; 
(e) k±i = —0.3, k±2 = 0.3; (f) k±i = —0.5, k±2 = 0.5. Maxima of figures correspond to probability 
densities: (a) 9.1 x 10"^^ (b) 6.5 x 10"^^ (c) 8.36 x 10"^^ (d) 8.4 x 10"^^ (e) 9.1 x 10"^^ (f) 
2.8 X 10-20. 
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FIG. 4: Role of electron-core interaction. Correlated spectra of Ar (linear scale) at 7x 10^^ W/cm^, 
800 nm are calculated using Eq. (19) with ryl = ''o = 0.985 a.u. (7 = 1.373). Both electron- 
elecron and electron-core interactions are included. Spectra are shown for different values of A;j_ (in 
a.u.) for both electrons: (a) /c^i = k^2 = 0; (b) = 0, k±2 = 0.2; (c) k±i = —0.1, k±2 = 0.2; (d) 
k±i = —0.2,k±2 = 0.3; (e) k±i = — 0.3,A;_l2 = 0.3; (f) k±i = —0.5,k±2 = 0.5. Maxima of figures 
correspond to probability densities: (a) 2.8 x 10"^ (b) 3.3 x 10"^ (c) 6.1 x 10"^ (d) 1.5 x 10"^, 
(e) 4.4 X 10-^ (f) 7.7 x 10"^ 
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